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Abstract 

Matrix elements of spin-dependent twist-four operators are extracted from 
recent data on the spin-dependent g\ structure function of the proton and 
deuteron in the resonance region. We emphasize the need to include the 
elastic contributions to the first moments of the structure functions at Q 2 < 2 
GeV 2 . The coefficients of the 1/Q 2 corrections to the Ellis- Jaffe sum rules 
are found to be 0.04 ± 0.02 and 0.03 ± 0.04 GeV 2 for the proton and neutron, 
respectively. 
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There has been much activity in recent years surrounding various deep-inelastic spin sum 
rules, such as the Bjorken |l] and Ellis- Jaffe sum rules 0], which test our understanding of 
the spin structure of the nucleon, as well as our ability to calculate higher-order radiative 
corrections in QCD. An intriguing issue in the study of these sum rules at moderate values 
of the momentum transfer squared, Q 2 (say Q 2 ~ 0.5 to 3 GeV 2 ), is that of higher-twist 
corrections. In this paper, we extract the twist-four matrix elements from recent data taken 
by the E143 Collaboration at SLAC |J on the proton and deuteron g\ structure functions, 
and compare these with some recent theoretical estimates. 

Higher-twist corrections to the polarized deep-inelastic sum rules were first studied by 
Shuryak and Vainshtein [|]]. The coefficient functions were recalculated in Ref. [||] and 
confirmed in Refs. || and [[?]]. There is a long list of calculations and estimates of the non- 
perturbative higher- twist matrix elements in the literature P,p|,|3|-|i"2|| . The interplay between 
the perturbation series at high orders and higher-twist matrix elements was first discussed by 
Mueller in Ref. ||13|| . The question concerns the precise definition of higher- twist corrections, 
as these are related to the procedure of how to regularize an asymptotic perturbation series. 
A concrete proposal of separating the perturbative and non-perturbative contributions was 
suggested in Ref. fb3fl . 

Very recently, the E143 Collaboration |3j published the first data on the first moments 
of the gi structure functions of the proton and neutron (the latter being extracted from 
deuterium data) in the resonance region (at Q 2 = 0.5 and 1.2 GeV 2 ). According to the 
phenomenon of parton-hadron duality, these data have a direct bearing on the size of the 
higher-twist matrix elements. From the successful phenomenology of the QCD sum rule 
method [HJ, power corrections in an operator product expansion (OPE) have direct control 
over the structure of the low- lying resonances. In fact, much work has appeared in the litera- 
ture on the calculation of low-energy hadronic observables in terms of vacuum condensates. 
Of course, the parton-hadron duality also allows determination of the same condensates 
from the resonance masses and widths, if these are known to sufficient accuracy. In our 
case, we shall employ the latter approach to duality: namely, resonances (including the nu- 
cleon elastic contribution) fix the higher-twist matrix elements. In the case of unpolarized 
deep- inelastic scattering, a similar analysis was performed in Ref. [[16] . Here we shall use the 
newly measured data from Ref. to extract the spin-dependent twist-four matrix elements. 

The first moment of the g\ structure function of the nucleon is defined as: 

T N {Q 2 )= ( l dxg»{x,Q 2 ) , (1) 
Jo 

where N = p or n, and the upper limit includes the nucleon elastic contribution. The 
inclusion of the elastic component is critical if one wishes to use the OPE to study the 
evolution of the sum rule in the moderate Q 2 region [I7|. Note that the application of the 



OPE requires a product of two currents, which arises in deep-inelastic scattering only after 
summing over all final (including elastic) hadronic states. For unpolarized scattering, this 
point was emphasized some time ago by J. Ellis [18]. In Ref. f[9 |, attempts were made to 



extract the twist-four matrix elements from the first moment of gi, without inclusion of the 
elastic contribution. In light of the logic behind the OPE, this procedure is clearly incorrect. 
Furthermore, the OPE is known to break down at low Q 2 (Q 2 ^ 0.5 GeV 2 ); hence the 
Drell-Hearn-Gerasimov sum rule |20| at Q 2 = has no obvious bearing on the size of the 
twist-four contributions IB. 
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In Figs.l and 2 we show the first moments T P (Q 2 ) and T n (Q 2 ), respectively, at three 
different Q 2 values (Q 2 ^ 3 GeV 2 ) (the data are from Refs. ||^,^|). The data, shown by the 
open circles, contain contributions from the measured x-regions covered in the experiments, 
together with theoretical extrapolations into the unmeasured small- X region, and also from 
elastic scattering at x — 1. Note that the neutron points are obtained from the deuteron and 
proton data assuming a small (~ 5%) D-state admixture in the deuteron. Compared with 
the proton data, the neutron results are much poorer, and could be significantly improved 
in future measurements at Jefferson Lab |24| . One should note that the contributions from 
the x- region covered by the recent E143 analysis [Q] are not entirely model independent. 
A number of assumptions have been made regarding the contributions of the #2 structure 
function, the A resonance, etc. |§. For the present analysis, however, we simply adopt the 
data as published in Ref. ||, keeping in mind that future experiments may help to clarify 
some of these assumptions. 

For the small-x extrapolation, we use a recent parameterization from Ref. [^] fitted 
to the global polarized deep-inelastic scattering data, which respects the Bjorken sum rule 
(gA = 1.257 and a^ LO (M|) = 0.117 ± 0.005). The parameterization has been constructed 
at a scale Qq = 1 GeV 2 , and evolution to different Q 2 values has been done including the 
complete NLO corrections [275|. At such small Q 2 values one might question the role of 



higher-twist contributions in any small-x extrapolation. However, previous experience from 
unpolarized deep-inelastic data tells us that higher-twist effects at small x tend to be rather 
small p7j . It is difficult to assign an error to the theoretical small-x extrapolation, as this 



is still a somewhat controversial, but interesting, subject which is currently under active 
study p8|-|32|. The data shown by the squares in Figs.l and 2 include only the inelastic 



contributions, as discussed above. 

The dotted curves in Figs.l and 2 represent the elastic contributions to the T N moments 
calculated in terms of nucleon form factors ||: 

r^(Q 2 ) = \f*{q 2 ) (F? (Q 2 ) + F 2 N (Q 2 )) , (2) 



where we have used the parameterization of Ff 2 (Q 2 ) from the fit of Ref. For the proton, 
the elastic contribution is negligible at Q 2 > 3 GeV 2 , and is about 10% of the inelastic at 
Q 2 = 2 GeV 2 . At Q 2 = 1 GeV 2 it is as important as the inelastic component, and below 0.5 
GeV 2 the elastic contribution becomes dominant. For the neutron, the elastic contribution 
peaks around Q 2 = 0.5 GeV 2 , and becomes quite small above ~ 1.5 GeV 2 and below « 0.1 
GeV 2 . 

According to the OPE, in the limit of large Q 2 ^> Aq CD , T n (Q 2 ) can be calculated via 
the twist expansion: 

r*W) = E ^§r, (3) 

r=2,4,--- *V 

where ^ is related to nucleon matrix elements of operators of twist < r. 

From the total T N (Q 2 ) one can obtain the higher-twist component by subtracting the 
twist-two contribution, fi 2 , which can be written as a series expansion in a s : fi 2 (Q 2 ) = 
J2 n &n a s (Q 2 ). It is suspected that the coefficients C£ grow like n\ as n — > oo, and therefore, 



strictly speaking, fi 2 is n °t a well-defined quantity [T3|]. The uncertainty in regularizing the 
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divergent series is closely related to the precise definition of the higher-twist contributions. 
In this paper we maximally utilize the available perturbative calculations up to 0{a z s ) ||33|| , 
and define n 2 up to this order as the entire twist-two contribution. We will return to this 
point later. For three quark flavors, the three-loop result for the twist-two component of the 
proton and neutron first moments is given by p3f : 
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(4) 



where the ± refers to p or n. The leading-twist component in (^) is given in terms of 
the triplet and octet axial charges, qa and a 8 , respectively, and the quantity defined 
as the renormalization group invariant nucleon matrix element of the singlet axial current 
33| , 



£(Q 2 



IS 



oo). With the above value of gA, the leading-twist contribution 
calculated with the values a 8 = 0.579 ±0.025 [31], extracted from weak hyperon decays, and 
m 0.15 ± 0.12 obtained from the global fit of Ref. p5|. The results for the proton and 



neutron are shown in Figs.l and 2 by the solid curves, where the band range reflects the 
combined error in Soo and a,. For the latter we have taken a? LO (l GeV 2 ) = 0.45±0.05 M 



Subtracting from the open circles in Figs.l and 2 the leading- twist component, one obtains 
the pure higher-twist contribution, denoted by the full circles. Note that the increased error 
bars on the higher-twist points simply reflect the uncertainty in fi 2 {Q 2 )- Thus the higher- 
twist contribution and the uncertainty in a s are correlated at intermediate values of Q 2 
(0.5 < Q 2 < 2 GeV 2 ). The higher-twist contribution can thus be reliably extracted only 
when a s is more accurately determined from other sources. 

Turning now to the higher-twist contributions to T N (Q 2 ), the 1/Q 2 term in Eq.(||) is a 
sum of three terms: 



^(Q 2 ) = \M 2 (a N 2 {Q 2 ) + Ad N 2 {Q 2 ) - 4f 2 N (Q 2 



(5) 



The a 2 component, being given by the second moment of the twist-two part of the polarized 
gf structure function, 



dx x g x (x, Q ) 



(6) 



arises from the target mass correction || , and can be evaluated straightforwardly from global 
parameterizations of the leading-twist part of g±(x,Q 2 ). The twist-3 correction in Eq.(^) 
can be extracted from the leading-twist contributions to the following moment of the g^ 
and g 2 structure functions: 



dx x 



(2^(x,Q 2 



(7) 



Recently the d 2 coefficients of the proton and deuteron have been determined in Ref. B5 
at an average value of Q 2 = 5 GeV 2 , with the results: d 2 = 0.0054 ± 0.0050, and d 2 = 
0.0039 ±0.0092. To a good approximation the neutron d 2 can be obtained from the relation: 
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d% = 2d>2 / (1—3/2 uj£>)— d%, where ujd is the deuteron D-state probability (more sophisticated 
treatments |36]| which account for binding and Fermi motion effects are small in comparison 



with the present error bars). With u D = 5% this gives d% = 0.0030 ±0.020 at Q 2 = 5 GeV 2 . 



Without any D-state correction |37j] the value would be around 25% smaller. To obtain d 2 



at a different Q 2 one can use the leading logarithmic evolution as computed in Refs. p|,|38 
In principle, the anomalous dimensions for the singlet and non-singlet components differ, so 
that one would require separate knowledge of the singlet and non-singlet matrix elements 
to perform the evolution. In practice, however, for 3 flavors the one-loop singlet and non- 
singlet anomalous dimensions turn out to be very similar |3]J3£|, and at the present level of 
accuracy these differences can be safely ignored. 

Lastly, the twist-four contribution, f^, to /i^ is defined by the matrix element: 

2f 2 N (Q 2 )M 2 S» = $> 2 (P, S \g^ f F^^ f \ P, S) , (8) 
/ 

where S M is the nucleon spin vector, and the gluon field-strength tensor is F^ v = 
1/2 e^^Fap, with e 0123 = +1, and g appears in the covariant derivative as _D M = <9 M + igA fl . 

In the remainder of this paper, our focus will be on the extraction of this matrix element 
from the data in Figs.l and 2. To achieve this, one must first subtract the leading-twist, 
target mass, and the twist-3 contributions from the T N (Q 2 ) data, and define: 

AI^(Q 2 ) = T^(Q 2 ) - ^(Q 2 ) - 1^ (a»(Q 2 ) + <(Q 2 )) . (9) 

In Figs. 3 and 4 the extracted AT N (Q 2 ) values for the proton and neutron are shown, respec- 
tively, as a function of 1/Q 2 - Theoretically, the data in Figs. 3 and 4 represent contributions 
from f^(Q 2 ) as well as from r = 6 and higher twists: 

AT N (Q 2 ) = -\^{Q 2 ) + E ^ ■ (10) 

The twist expansion is believed, however, to be controlled by a scale related to the average 
transverse momentum of quarks in the nucleon [39], typically of the order 0.4-0.5 GeV [|5,39 



Therefore one can reasonably expect that the role of r > 6 effects should not be significant 
for Q 2 > 1 GeV 2 , and not overwhelming for Q 2 > 0.5 GeV 2 . 

Finally, the matrix elements f% can be extracted from the data points in Figs. 3 and 
4 at Q 2 > 1 GeV 2 by neglecting the higher- twist terms in Eq. (|T0|) . The Q 2 evolution of 
f^(Q 2 ) is also taken into account at leading logarithmic order [^jp8|. This logarithmic Q 2 
dependence results in the slight deviations in the curves in Figs. 3 and 4 from linearity. At 
a scale of Q 2 = 1 GeV 2 we find: 

/f = -0.10 ±0.05 , (11) 
/£ = -0.07 ±0.08 . (12) 

To determine the effects of the higher-order terms in the coefficient functions of the twist- 
two contributions, we have also performed the extraction by neglecting the 0(a^) terms in 
Eq. (f|) . The central value of the twist-four matrix element for the proton is then reduced from 
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—0.10 to —0.08, while for the neutron it increases slightly, from —0.07 to —0.08. Including 
the a% and contributions, one can finally determine the 1/Q 2 correction to the Ellis- Jaffe 
sum rules: 

\i\ = (0.04 ± 0.02) GeV 2 , (13) 
H2 = (0.03 ± 0.04) GeV 2 , (14) 

at the scale Q 2 = 1 GeV 2 . 

The central values in Eqs.([Tl|) and (|l^) seem to suggest that the isoscalar combination 
of the twist-four matrix elements is much larger than the isovector combination. One might 
suspect therefore that the singlet twist-two contribution to the sum rule obtained from the 



global fit p5| is too small. To investigate the effect that a larger value of Eqo would have on 
the twist-four matrix elements, we have reanalyzed the data using T,^ « 0.3 as the central 
value [P2"] . The effect is a reduction of /f to ~ —0.05, and f£ to ~ 0.0, which would then 
lead to similar values for both the isotriplet and isosinglet combinations. Of course these 
values are still consistent with the results in Eqs.(|Tl"D and (|T2]) within the errors. In principle, 
one could eliminate the dependence on Soo by considering only the isovector combination 
r p — r n , thereby reducing significantly the uncertainty in the isovector twist-four matrix 
element. Unfortunately, the error associated with the neutron data is largely experimental, 
so that the final proton — neutron moment would have an error which is as large as that for 
the neutron points in Fig. 4. 

The values determined in Eqs.flTTD and (|12|) can be compared with several model calcula- 
tions of the twist-four matrix elements in the literature. The first estimates of were made 
using QCD sum rules. The result from Ref. @ is /f = 0.050±0.034 and / 2 n = -0.018±0.017, 



while that from Ref. P] is /f = 0.037±0.006 and / 2 n = 0.013±0.006. Alternative estimates 



of the t = 4 matrix elements were made using the MIT bag model Hf40|. The result there, 
evolved from the bag scale up to Q 2 ~ 1 GeV 2 , was found to be /f = —0.028 and = 0. 

The results obtained in this work will be improved as more experimental information 
becomes available in future. The error on the data points in Figs. 3 and 4 come mainly 
from the uncertainty associated with the value of a s , and the singlet axial charge Eqo, when 
subtracting the twist-two contribution from the total T N (Q 2 ). Therefore better knowledge of 
the twist-two part of the structure function at higher Q 2 , and a more accurate determination 
of ot s from other experiments, will be valuable in pinning down the higher-twists at low 
Q 2 . Certainly more data points are needed in order to establish a clearer trend of the Q 2 
dependence at moderate values of Q 2 (Q 2 ~ 0.5 — 3.0 GeV 2 ). In particular, the neutron 
data points are irregular and should be confirmed in subsequent experiments. In this respect, 
future experiments at Jefferson Lab and other facilities can contribute much to our present 
understanding of the twist-four matrix elements of the nucleon. 
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FIG. 1. Q 2 dependence of T P (Q 2 ). The solid curves represent the twist-two part of T P (Q^ 



while the dotted is the elastic component, as parameterized in Ref. [21|. The squares denote the 
inelastic contribution extracted from the E143 experiment |J, the open circles include also the 
elastic piece, while the full circles are the pure higher-twist contributions. 
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FIG. 2. Same as in Fig.l, but for the neutron. 
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FIG. 3. The 1/Q 2 dependence of the higher-twist contribution AT P (Q 2 ) defined in Eqs.(^) and 
(|lT)|). The curves correspond to f% = —0.10 ±0.05 at Q 2 = 1 GeV 2 (the solid represents the central 
value, while the dotted curves indicate the error range on f%). 
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FIG. 4. Same as in Fig. 3 but for the neutron. The curves correspond to f£ 
Q 2 = 1 GeV 2 . 
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